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Advanced techniques for manipulation of internal states, standard in atomic physics, are demon- 
strated for a charged molecular species for the first time. We address individual hyperfine states of 
ro- vibrational levels of a diatomic ion by optical excitation of individual hyperfine transitions, and 
achieve controlled transfer of population into a selected hyperfine state. We use molecular hydrogen 
ions (HD^) as a model system and employ a novel frequency- comb-based, continuous- wave 5 pm 
>! laser spectrometer. The achieved spectral resolution is the highest obtained so far in the optical 

^ ^ domain on a molecular ion species. As a consequence, we are also able to perform the most precise 

^0-^ test yet of the ab-initio theory of a molecule. 

<N 

H Cold trapped molecules [TJ |2] currently represent an intense field of activity relying on sophisticated methods of 

^^ molecule production, translational and internal cooling, spectroscopy and sensitive detection. Many applications, such 

^^ as chemical reaction studies fSjIJ, tests of molecular quantum theory [5j, fundamental physics [6, 7J and quantum 

-^ computing [8j would benefit strongly from the availability of advanced manipulation techniques, already standard in 

^Nl atomic physics. These are not straightforward for molecules, and for charged molecules have not been demonstrated 

yet. Production methods for molecular ions (usually by electron impact ionization) and, if heteronuclear, their 

interaction with the black-body radiation of the surrounding vacuum chamber, usually lead to significant population 

of a substantial number of internal states. A first, important step in the manipulation of internal states of molecular 

ions is population transfer between rotational states (heteronuclear molecules usually being cold vibrationally, i.e. 

C are all in the v = ground vibrational state). It has been demonstrated that a significant fraction (ca. 75%) of an 

Q ensemble of diatomic molecular ions can be pumped into the vibrational and rotational ground level (v = 0, TV = 0) 

^lOj, see Fig. [I] 

For a general diatomic molecule, however, this pumping is usually not capable of preparing molecules in a single 



Oh 



C^ 



^ quantum state, because spin interactions generate a hyperfine structure with several states in each ro- vibrational level. 

. ^-H For example, a diatomic molecule with one unpaired electron (sg = 1/2), and nuclei with nuclear spins /i = 1/2, /2 = 1 

^^ (such as HD+) has 4 hyperfine states in zero magnetic field if the rotational angular momentum A/" = 0, but 10 if 

Cj A^ = 1, and 12 if A" > 2, see Fig. [2] a. The ability to address selectively molecules in one particular hyperfine state 

O (or even in a single quantum state with a particular magnetic quantum number J^) and to transfer molecules from 

•—I one hyperfine state to another are clearly important tools of a molecular quantum toolbox that can be part of a full 



<N 

> 






quantum state preparation procedure. 

Complicating the addressing, the number of strong transitions between two given ro-vibrational levels 
('U, A^), {v' ^ N') is equal to the larger of the two numbers of hyperfine states, i.e. potentially high, and with only small 



i/^ differences in transition frequency. Fig. [2| b shows as an example the case of the fundamental vibrational transition 
T-H ('U = 0, A" = 0) ^ {v' = l^ N' = 1) in HD^, where 10 strong transitions occur over a range of about 60 MHz (pT|[T2]). 

CN Addressing a single hyperfine state in a multi-spin molecule thus requires a spectroscopy that can resolve individual 

^^ "hyperfine" lines in the spectrum. 

^^ In this work, our approach is based on one-photon laser excitation of the fundamental vibrational transition (0, 0) -^ 

CN (1, 1) at the wavelength A/, see Fig. fl] The relatively low transition frequency. A/ > 2.5 /im for diatomics, in 

^"^ combination with the low secular kinetic energy ksTsec achievable by sympathetic cooling, yields a Doppler broadening 

^ Aujj of the transitions that is smaller than many line spacings. This provides the desired quantum state selectivity for 

addressing some of the hyperfine states, using strong transitions. For our test case HD+, A/ = 5.1 /im, Tsec — 10 mK, 
Aud :^ 3 MHz. Additionally, excitation of weak transitions (which violate the approximate selection rules AF = 0, 
AS = 0, see Fig. [2]b), provides selectivity for all hyperfine states, since for these transitions the frequency spacings 
are larger. Compared to the use of a pure rotational excitation (0, 0) -^ (0, 1) or a microwave transition within a 
ro-vibrational level, the use of a vibrational transition has the advantage that the excitation may be followed by a 
much faster spontaneous decay, either back into the ground ro-vibrational level (rate approx. 6 s~^, here) or into the 
(relatively long-lived) rotational level {v" = 0, A"'' = 2) (rate approx. 12 s~^, here). This allows reasonably rapid 
pumping of the molecule (possibly after repeated absorption and spontaneous emission events) into another long-lived 
state, a necessary condition for efficient quantum state preparation, as shown below. 

The experiment is performed on ensembles of HD+ ions trapped in a linear quadrupole radio- frequency trap 
(14.2 MHz), sympathetically translationally cooled by co-trapped, laser-cooled Beryllium atomic ions [13j and rota- 
tionally cooled by lasers. Our laser system consists of four subsystems: the A/ = 5.1 pm laser spectrometer referenced 
to an atomic frequency standard [14j (see Supplementary materials), a reliable, frequency-stabilized fiber-laser-based 



313 nm laser for cooling of Beryllium ions ^5], a rotational cooling laser system (a Xp = 5.4 pm quantum cascade laser 
and, for part of the measurements, a A^ = 2.7 pm diode laser), and a pair of lasers (A, A') for resonance-enhanced 
multi-photon dissociation (REMPD), see Fig. IT] Rotational cooling |9j is a crucial tool here, as it significantly in- 
creases the fractional population of molecules in the lower ro- vibrational level ('U = 0, A/" = 0), from ca. 10% to 60 - 
75%. The possibility provided by our laser system to measure the HD+ fundamental vibrational transition frequency 
/ = c/Xf also allows us to perform a precise comparison with the ab-initio theory of the molecular hydrogen ion. 

Hyperfine-resolved ro- vibrational transitions are induced by the Xf = 5.1 pm laser tuned precisely to individual 
transitions. In order to show that this is actually achieved, we detect (without hyperfine state selectivity) the 
population of the goal vibrational level {v^ = 1, A^' = 1), by 1+1' resonance-enhanced multi-photon dissociation 
(REMPD)[16], see Fig. [l] Together, this represents a three-photon (l+l'+l") REMPD process. The reduction of the 
HD+ number by the REMPD process is our spectroscopic signal [16j. 

Figure [3] shows the obtained hyperfine-state resolved spectrum. All theoretically predicted and addressed hyperfine 
transitions were observed; only the (nominally strong) transition S4 is barely detected, for unknown reason. Transition 
W4, which originates from a lower hyperfine state with only a small fractional population (J = 0, thus statistically 
containing only 5 - 6% of all molecules), could only be made clearly evident using a preceding hyperfine pumping 
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Figure 1: Schematic energy level scheme of HD^ with transitions relevant to this work. Hyperfine structure is shown 
schematically only for the {v — 0, N — 0) and (1, 1) levels as lines, but is implied for all other levels as well (thick bars). Full 
lines: laser-induced transitions, dashed lines: some relevant spontaneous emission transitions; dotted lines: some black-body 
induced transitions. The spectrally narrow wave A/ selectively excites molecules from a particular hyperfine state {v = 0,N = 
0,F, S, J) to a single hyperfine state (1, 1, F\ S\ J'). Quantum state preparation is performed by irradiating alternatingly 
the appropriately tuned wave A/ and Xp , in conjunction with spontaneous emission from the level (1, 1). Resonant laser 
radiation at A^ and nonresonant radiaton at A^^ is used to detect that hyperfine-state-selective excitation to the (1,1) level has 
occurred, by transferring the excited molecules to the electronically excited molecular state 2pcr from which they dissociate. 
Initially, rotational cooling is performed by radiation at Ap and A^ . The level energy differences are not to scale. The waves 
at A^, A^^, Ap, Ap have large spectral linewidths and do not excite hyperfine state-selectively. 
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Figure 2: a (left): Energy diagram of the hyperfine states and main electric-dipole transitions in zero magnetic field, b (right): 
Stick spectrum of the transitions (in zero field; values of the squared transition moment (f are normalized to the strongest 
transition). The states are labeled by the quantum numbers {F, S, J). Weak transitions are shown in pink [12 . Very weak 
transitions are not shown, except for P4 (dashed). The "spin-independent" transition frequency fo^th is the value if nuclear and 
electron spin were zero. SO, SI, S2, S3, S4, S5, S6, Wl, W2, W3, W4, P2 are transitions studied here ("W, S, P" mean "weak", 
"strong", "pumping", respectively). All were observed except S4. Wl, W3 are the transitions used here to achieve population 
transfer from the hyperfine states (v = 0, N = 0, F = 1, 5* = 2, J = 2, Jz) (empty circle) and (0, 0, 1, 1, 1, J^) into the 
hyperfine state (0, 0, 1, 0, 0, J'J) (filled circle). PI, P2, P3, P4 are proposed optical pumping transitions (with indicated 
polarizations) for preparation of the molecule in the single quantum state (0, 0, 1, 2, 2, Jz = +2) (one of the Zeeman states in 
the open circled hyperfine state). 



step (see below). Each of the 4 hyperfine states of the lower level was selectively addressed, and 5 of the 12 upper 
level hyperfine states were selectively populated. We also observed the line at —10.2 MHz, which contains two nearly 
coinciding transitions S2, S3, but originating from diflFerent ground hyperfine states. The remaining strong transitions 
(including the line marked PI in Fig. [2| were also observed, but their small spacing prevents complete resolution, 
and they are not reported in the Figure]3] 

We demonstrate hyperfine state manipulation by optical pumping of individual hyperfine state populations into a 
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Figure 3: Observed hyperfine spectrum of the (^7 = 0, N = 0) -^ (1, 1) fundamental ro- vibrational transition in cold trapped 
HD^ ions. The effective intensity times irradiation duration product of the 5.1 pm radiation varied from line to line, and 
was adapted to avoid saturation. Brown lines are the result of fitting fo,exp and the individual line amplitudes for 9.5 mK 
temperature, and an average magnetic field of 0.8 G. The sticks are for illustration purpose and show the theoretical squared 
transition dipole moments for the Zeeman components at 0.8 G, assuming exciting radiation polarized at 45 degrees to the 
magnetic field. They are scaled by different factors for presentation purpose. Color coding is as in Fig. [2] S4, S6 were taken 
at high intensity-irradiation time product, S5 at a lower value. The side peak of SI is probably due to an ion micromotion 
sideband of S2/S3. The W4 line required hyperfine state optical pumping for its detection (see Fig. p|. 



goal state (filled circle in Fig. [2|. As a goal state we choose (v = 0, A/" = 0, F = 1, S* = 0, J = 0, J;^ = 0) which is 
non-degenerate (J = 0) and thus a single quantum state. After rotational cooling, we apply the following sequence 
twice: Wl line (3 s), rotational re-pumping (Ap and A^ simultaneously for 5 s), W3 line (3 s). A final 10 s of rotational 
re-pumping is performed before the spectroscopic excitation. The Wl and W3 transitions excite population from 
two initial hyperfine states (0, 0, 1, 2,2, J^), (0, 0, 1, 1, 1, J^) (without Jz selectivity) into the same hyperfine state 
(1, 1, 1, 0, 1, J'J). This state has dominant spontaneous decay to the goal state (green line S3 in Fig.[2|. We find clear 
evidence that this hyperfine state preparation is taking place by observing the transition W4 starting from the goal 
state by REMPD, see Fig. [4j This transition is not observable in our experiment without the preparation procedure, 
since then the population in the lower hyperfine state is too low. 

Our hyperfine-state resolved spectrum represents the highest-resolution optical spectrum of any molecular ion 
so far [51 [171 HH]- This enables an accurate comparison of experimental frequencies with ab-initio theory. 
We obtain two hyperfine state separations in the ground state, A/o,o,a = (^(0, 0, 1, 2, 2)— £^(0, 0, 1, 1, l))//i and 
^/o,o,6 = (^(0? O7 17 17 1)~^(07 O7 17 O7 0))/^7 from the measured transition frequency combinations f{Sl) — /(W2), 
f{wi)-f{Wl), f{P2)-f{S0), and from f{W4)-f{Sl), /(P3) -/(P2), respectively. In addition, two excited state 
splittings, A/1,1,, = (F(l,l,l,2,l)-F(l,l,l,l,l))//i and A/i,i,rf = (F(l, 1, 1, 1, 1)-F(1, 1, 1,0, l))//i are similarly 
obtainable by suitable frequency combinations. A fit of these hyperfine state separations to the data (fitting also 
the spin-independent frequency) yields agreement with the ab-initio results (A/0,0, a^ ^/o,o,67 ^/i,i,c7 ^fi,i,d)th — 
(130.60(1), 82.83(1), 113.33(1), 71.68(1)) MHz [11,119], with deviations {exp. - theory) of (-0.22(0.13), 0.28(0.38), 
—0.13(0.15), 0.27(0.19)) MHz. The two measured hyperfine separations of the ground state also allow determining the 
two hyperfine constants ^4(0, 0), ^5(0, 0) fTT] which fully describe the hyperfine structure of the ground state [12J. Our 
fit yields (^4(0, 0), ^5(0, 0)) = (906(17), 142.33(25)) MHz, whereas the theory values are (925.38(1), 142.29(1)) MHz 
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Figure 4: Demonstration of hyperfine state manipulation. The transition W4 shown here is observed only when hyperfine 
optical pumping is implemented. This transition represents the excitation from a single quantum state, (0, 0, 1, 0, 0, Jz = 0). 
Data shown was taken alternating measurements preceded by hyperfine optical pumping (upper data points joined by line) 
and not (lower data points). The intensity of the 5.1 pm laser was set to its maximum both during hyperfine pumping on 
the Wl and W3 transitions and subsequent detection of the W4 transition. Irradiation time on the W4 transition was 3 s. 
Rotational cooling by the 2.7 pm and 5.5 pm laser was used. The zero level corresponds to the relative decrease measured 
when the 5.1 pm spectroscopy laser was blocked. The three sticks show, for illustration purposes, the theoretical transition 
frequencies and strengths in a 0.8 G magnetic field and radiation polarized at 45 degree to the magnetic field. The shift of the 
central component is —0.05 MHz relative to the zero-field frequency. 



Assuming instead that the hyperfine energies are given by the theoretical values (this assumption being strenghtened 
by the agreement of hyperfine theory and experiment for large-v levels pT| |20j), we can fit an overall frequency 
correction to the spectra Wl, W2, W3, SO, SI, S2+S3, W4, P2. We obtain the spin-independent frequency fo,exp = 
58 605 052.00 MHz, with combined statistical and systematic error of 0.064 MHz (see Supplemental Materials). 
The theoretical value is fo^th = 58 605 052.139(11) (21) MHz, where the first error is due to the uncertainty of the 
fundamental constants and the second is the theoretical error in the evaluation of the QED contributions p!9l [2T | [22] . 
The difference between exprimental and theoretical resuls is —2.0 times the combined theoretical plus experimental 
error. The relative experimental uncertainty of 1.1 x 10~^ represents the most accurate test of molecular theory to date. 
In particular, our measurement is the first molecular measurement sufficiently accurate to be explicitly sensitive to 
the QED contributions of order a^ (relative to the nonrelativistic contribution to the transition frequency), calculated 
as 0.109(21) MHz for the transition studied here. 

Based on the technique demonstrated here, we can propose a realistic optical pumping procedure for preparing 
most of the population in a single quantum state (v, N, F, S, J, J^), i.e. with well-defined projection of the total 
angular momentum. Under typical conditions, the relative statistical occupation of any individual quantum state in 
{v = 0^ N = 0) is only ~ (1/12) x (60% — 75%) ^ 5% — 6% under rotational cooling by a single laser (Ap) or two lasers 
(Ap, Ap. Exciting sequentially the four transitions P4 [(0, 0, 0, 1, 1) -^ (1, 1, 1, 2, 2)], P3, P2, P1[(0, 0, 1, 2, 2) -^ 
(1, 1, 1, 2, 3)] in a weak magnetic field and with polarizations chosen as indicated in Fig.|2]a will cause transfer of the 
population of all Zeeman quantum states of (0, 0) to the single Zeeman quantum state (0, 0, 1, 2, 2, J^ = +2), via 
spontaneous emission processes from (1, 1, 1, 2, J' = {1, 2}, J^), which dominant ly occur on strong transitions (red 
lines in Fig. [21 a). These excitations should be interleaved with rotational cooling (lasers Ap, A^, which also serves as 
repumper following spontaneous decay into {v" = 0, N^' = 2). The optical pumping procedure should take a few ten 
seconds and lead to 60% - 70% fractional population in the goal state. 

In summary, we have shown that it is possible to address and prepare individual hyperfine states in cold, trapped 
diatomic molecular ions even in presence of a complex spin structure. A mid-infrared laser spectrometer controlled by 
an atomic standard-referenced frequency comb, and sufficiently low ion kinetic energies were two important require- 



ments. The observed, Doppler-limited, transition linewidths (3 MHz) are the lowest obtained to date on a molecular 
ion species in the optical domain (note that they scale as (molecule mass)~^/^). We also observed, for the first time 
to our knowledge, weakly allowed hyperfine transitions using optical excitation. As one application, we were able to 
directly determine the population fraction of molecules in particular hyperfine states. The largest value we found was 
19%, clearly indicating the effectiveness of our rotational cooling. We also demonstrated excitation of a transition 
from a single quantum state. Since our test molecule HD+ is the simplest heteronuclear molecule and is excited 
from the ro- vibrational ground state, this study represents the first precision measurement of the most fundamental 
electric-dipole allowed ro- vibrational transition of any molecule [17J. A comparison of theory with experiment showed 
that (i) the hyperfine energies of small- v, N ro- vibrational levels agree within deviations of less than 0.3 MHz and (ii) 
the spin- independent energy agrees within 2 times the error of 1.1 x 10~^. 
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Supplemental Materials 
1. General features of hyperfine structure 

We consider a diatomic molecule with total electron spin Sg, nuclei with nuclear spins /i, /2, in a rotational level 
N. The spin degeneracy is lifted (only partially in zero magnetic field) into a number of hyperfine states by the 
electron spin - nuclear spin interaction (~ s • Ii) and/or the electron spin - rotation interaction (~ s • N). In HD+, 
the (approximate) quantum numbers (F, S^ J) correspond to the couplings F = s + Ip, S = F + Id, J = S + N. The 
electron spin - nuclear spin interactions determine the main splittings in a given ro-vibrational level, resulting in 4 
singlets (if TV = 0) or in 4 multiplets (if N y^ 0). In the latter case, the electron-spin rotation interaction determines 
the splitting strength within the multiplets pTj. 

The line splittings in the transition spectrum, Fig.[2]b, arise because in the lower and upper levels (i) the strengths of 
the electron spin-nuclear spin interactions and of the electron spin - rotation interaction differ, and (ii) the rotational 
angular momenta N^ N' differ. The density of lines in the spectrum is essentially independent of whether the transition 
is a fundamental vibrational transition, an overtone transition (A^' = v^ — v > 1) (e.g. [E]), or a pure rotational 
transition {Av = 0). 

2. 5.1 pm laser spectrometer 

The spectroscopy radiation was generated by frequency-mixing in a nonlinear optical crystal two near-infrared lasers 
which are individually referenced to a conventional frequency comb [14J. The generated wavelength, 5.1 pm, extends 
the coverge provided so far by sources with absolute frequency measurement capability [1, 2j. 

We use two continuous- wave lasers, a NdiYAG laser (1064 nm, ca. 6 W) stabilized to a Doppler-free resonance in 
molecular iodine via its second-harmonic wave at 532 nm, and a 1344 nm home-built external cavity quantum-dot 
diode laser (QD-ECDL). The main part of the 1344 nm wave is resonantly enhanced in a ring resonator containing a 
periodically poled MgOiLiNbOs crystal with poling period 25.4 pm, appropriate for difference frequency generation 
of 5.1 pm radiation. The 1064 nm wave is focused and overlapped with the 1344 nm wave in the crystal, but not 
resonated. Although the crystal strongly absorbs the generated 5.1 pm radiation, an output power of up to 0.1 mW 
is generated. Small fractions of the 1064 and 1344 nm waves are sent via an unstabilized single-mode optical fiber 
to another laboratory containing the femtosecond Tiisapphire frequency comb, where both frequencies are measured. 
The 1344 nm laser is actively frequency- stabilized to a mode of the frequency comb. The 1064 nm wave frequency is 
continuously measured by the comb. The comb itself is stabilized to a hydrogen maser or ultra- stable quartz oscillator, 
both referenced to GPS. 

Tunable, frequency-stable 5.1 pm radiation is produced by amplitude modulating the diode laser output wave with 
an integrated-optic modulator at a variable radio-frequency Q before sending the wave to the resonator. This generates 
two sidebands, and the resonator length is stabilized to resonate one of the two, generating a single frequency at 5.1 
pm. By changing the modulation frequency Q^ the 5.1 pm radiation can be smoothly and precisely frequency-tuned 
over 460 MHz. We determined the spectral width of the 5.1 pm radiation as Tf c^ 0.68 MHz. The frequency instability 
(Allan deviation) of this radiation is constant at ca. 23 kHz for integration times r up to 10 s, dropping to 4 kHz 
at r = 800 s. Both linewidth and instability are significantly lower than the Doppler width of the HD+ hyperfine 
transitions studied here. The absolute frequency uncertainty of the 5.1 pm radiation is less than 10 kHz when averaged 
over 1 h. 

3. Procedures and analysis 

Each data point in a spectrum measurement is acquired from one loading of HD+ ions into a laser-cooled Be+ ion 
cluster. Fig. |5] shows an example. 

From molecular dynamics (MD) simulation of ion cluster fiuorescence images [3j we deduce an upper limit for the 
secular temperature Tsec ^ 15 mK. After loading, the molecules are rotationally cooled by Xp o:^ 5.5 /im radiation 
(and, for some measurements, also Xp' =2.7 /im radiation) for 30 s, transferring about 60% (75%) of the molecules 
into the {v = 0^ N = 0) ro-vibrational ground level. This population is spread over the four hyperfine states. Then, 
the rotational pumping laser waves are blocked and the 5.1 /im spectroscopy radiation is applied along the trap axis 
for 3 - 30 s, simultaneously with a 1475 nm diode laser (ca. 1 mW) driving the {v^ = 1^ N^ = 1) ^ {v" = 5, N^^ = 2) 
overtone transition, and a 266 nm laser (ca. 25 mW) that dissociates the molecules in the latter state. 

The number of trapped HD+ ions is measured by excitation of their secular motion before the rotational cooling 
and after the (partial) photodissociation. Excitation of the secular motion heats up all ions and therefore changes 



Figure 5: CCD image of a typical Beryllium ion Coulomb cluster with embedded HD^ ions (darker region along the long axis), 
as used in the experiment. 

the Be+ fluorescence rate, which is detected by a photomultiplier. The relative decrease in the photomultiplier signal 
is the signal of interest, the excited molecule fraction. As the spectroscopy source linewidth is 0.7 MHz, we made 
frequency steps of typically 0.8 MHz when scanning through the lines. One data point could be recorded on average 
every 3 min. This limited the number of data points per frequency value that could be taken. We performed between 
4 and 25 measurements per frequency point. Typically, one hyperfine transition required one day of measurement. 
Background level measurements were performed in the same way, with the 5.1 pm laser blocked. Alignment of the 
5.1 iim laser with the ion cluster varied from day to day, which was compensated by adjusting the power. 

The typical standard deviation of the relative decay data taken at a given optical frequency (not of the mean) is 
0.02. This is due primarily to the low number of HD+ ions loaded into the ion cluster (a few hundred) and the low 
relative population in a particular hyperfine state (only a fraction c^ (0.6 — 0.75) pj, where pj = (2 J + 1)/12 in a 
particular hyperfine state J, for A^ = 0), which is further reduced during the REMPD phase (Ap, Ap ' lasers off) by 
the competing black-body-radiation driven excitation (^' = 0, A^ = 0) ^ (0, 1). The overall small numbers lead to 
significant statistical fluctuations in the number of ions actually prepared into a particular initial hyperfine state after 
optical pumping, and in the fraction excited and dissociated. As a consistency check, we found that the maximum 
observed dissociated fractions (at high intensity and on resonance) did not exceed pj. For example, for line SO, we 
observed (19.4±1.4)% ion number reduction, for lines S2+S3, (16.4±2.8)%. Note that the values shown in Fig. [3J are 
lower because there, lower intensities were used in order to avoid saturation broadening. 

The smallest observed full-width-half-maximum linewidths are ::^ 3.0 MHz. Assigning this to be due to Doppler 
broadening only, the value yields an upper limit of 15 mK for the secular temperature, consistent with the MD 
simulations. 

The observable spectrum may be modeled, including the effects of finite secular temperature, finite laser linewidth, 
finite population in the lower level, the effect of unresolved Zeeman splittings by the non-zero magnetic field in the 
ions' region [12j, and black-body-radiation induced excitation to (0, 1). In order to simplify the analysis, in the 
experiment we chose laser power and irradiation times such that significant saturation and concomitant broadening 
of the signals was avoided, i.e. the linear absorption regime was maintained. We furthermore take into account: (i) 
from independent measurements we have some knowledge about the magnetic field: an upper limit of ca. 1 G [4J; (ii) 
the laser linewidth does not contribute strongly to the total linewidth, and we may therefore use an effective Doppler 
temperature for modeling the lineshape; (iii) the overlap between the focused spectroscopy beam (ca. 0.5 mm waist) 
and the ion cloud was not constant over the time span covering all hyperfine line measurements; (iv) the effect of 
black-body-induced excitation (0, 0) -^ (0, 1) represents a reduction in maximum observable signal. Thus, we fit an 
effective product of intensity and irradiation time to each transition spectrum in order to reproduce the signal levels. 
We used a simplified model of the magnetic field inhomogeneity in the ions' volume. We find a good fit for an average 
magnetic field of 0.8 G and (9.5 ± 1) mK temperature. The statistical error for the spinless frequency fo^th is obtained 
from a Monte Carlo simulation as 60 kHz. 

4. Systematic errors 

Due to the slow data rate of this experiment, measurements of systematics were not possible. However, the well- 
developed theory of the HD+ molecule allows relying on theoretical results to estimate upper limits for various 
systematic errors. Those due to light shift, black-body shift [6j, electric quadrupole shift [7J, and Stark shift [8, 9j, 
are theoretically estimated to be less than 10 kHz. The only potential significant effect is the Zeeman effect. The 
sticks in Fig. |3J[4]show the splittings expected for the typical magnetic field value in our trap [12j. These unresolved 
splittings could lead to a shift of the centers of the lines but we find that the weighted mean frequency of the 
magnetic components of any line shifts weakly with magnetic field (e.g. less than ca. 10 kHz for line SO at 1 G). We 



modeled the influence of the imprecisely known magnetic field on the fitted spin-independent frequency and estimate 
an error of 13 kHz from this effect. In order to obtain the spin-independent frequency from our measurements, we 
use the theoretical prediction of the hyperfine energies. The influence of their theoretical errors onto the fitted spin- 
independent frequency is conservatively assumed to be 10 kHz. The error of the optical frequency measurements is 
less than 10 kHz. Similar considerations are applied to the fit of the hyperfine splittings. 

5. Theoretical transition frequencies 

The theoretical frequency of any particular hyperfine transition arises from three contributions. The main con- 
tribution is non-relativistic, 58 604 301.269 MHz, determined by the solution of the Schrodinger equation [21j and 
whose error (11 kHz) comes dominantly from the uncertainties of the electron-proton and proton-deuteron mass ratios 
TTie/mp^ and rrip/md [5^. A second contribution, 750.870 MHz, is from relativistic and QED effects [22j, with estimated 
theoretical uncertainty of 21 kHz [19j. The sum of these two contributions is the spin-independent frequency fo,th- 
The third is the hyperfine contribution [llj, e.g. -38.685 MHz for transition SO. Improved values for the hyperfine 
constants and Bethe logarithm communicated recently by V. Korobov have been used here. The influence of these 
improved values as compared to the published ones on the results presented here is small in comparison with the 
experimental uncertainties. The theoretical uncertainty of the hyperfine contributions to the transition frequencies is 
less than 10 kHz |19j . 

6. Hyperfine state preparation 

The re-pumping suggested acts predominantly on the strong, non-spin-state changing transitions. P2 is a weak 
transition, P3, P4 are very weak transitions. The transition dipole moments of the latter are so small [12j that a 
quantum cascade laser (typical output power level 10 mW) could be used for obtaining usefully large rates. Such a 
laser could be frequency-stabilized to a spectrometer of the type developed here, using the beat note with it as an 
error signal. 

An additional rotational cooling laser that depletes the (v = 0, L = 3) level by exciting from it to the (1, 2) or 
(2, 2) level could enable an even higher population fraction in the goal hyperfine state. 

In a similar way, any other Zeeman states of the ground vibrational level could in principle be populated to a high 
fraction. 

It is recognized that variations of the above scheme can also be used to achieve similar population fractions in 
individual quantum states of the {v = 0^ N = 1) level. 



[1] Wavelengths shorter than 4.4 jim were obtained in: P. Cancio, S. Bartalini, S. Borri, I. Galli, G. Gagliardi, G. Giusfredi, 

P. Maddaloni, P. Malara, D. Mazzotti and P. De Natale, Appl. Phys. B 102, 255 (2011) 
[2] For wavelengths near 10.6 ]im, see A. Amy-Klein, A. Goncharov, M. Guinet, C. Daussy, O. Lopez, A. Shelkovnikov, and 

C. Chardonnet, "Absolute frequency measurement of a SFe two-photon line by use of a femtosecond optical comb and 

sum-frequency generation", Opt. Lett. 30, 3320 (2005) 
[3] C. B. Zhang, D. Offenberg, B. Roth, M. A. Wilson, and S. Schiller, Phys. Rev. A 76, 012719 (2007) 
[4] J. Shen, A. Borodin, M. Hansen, and S. Schiller, subm. to Phys. Rev. A. (2011) 
[5] PJ. Mohr, B.N. Taylor, and D.B. Newell (2011), "The 2010 CODATA Recommended Values of the Fundamental Physical 

Constants" (Web Version 6.0). This database was developed by J. Baker, M. Douma, and S. Kotochigova. Available: 

http://physics.nist.gov/constants National Institute of Standards and Technology, Gaithersburg, MD 20899 
[6] J.-Ph. Karr, S. Kilic, and L. Hilico, J. Phys. B 38, 853 (2005); 
[7] A. K. Bhatia and R. J. Drachman, Phys. Rev. A. 61, 032503 (2000); 
[8] R. E. Moss and L. Valenzano, Mol. Phys. 100, 1527 (2002) 
[9] D. Bakalov and S. Schiller, to appear in Hyperfine Interactions (2012) 



